Singular point characterization in microscopic flows 
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We suggest an approach to microrheology based on optical traps in order to measure fluid fluxes 
around singular points of fluid flows. We experimentally demonstrate this technique, applying it 
to the characterization of controlled flows produced by a set of birefringent spheres spinning due 
to the transfer of light angular momentum. Unlike the previous techniques, this method is able to 
distinguish between a singular point in a complex flow and the absence of flow at all; furthermore 
it permits us to characterize the stability of the singular point. 
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The experimental characterization of fluid flows in 
micro-environments is important both from a fundamen- 
tal point of view and from an applied one, since for many 
applications it is required to assess the performance of 
microfluidic structures, such as lab-on-a-chip devices []]]. 
Carrying out this kind of measurements can be extremely 
challenging. In particular, due to the small size of these 
environments, wall effects can not be neglected Addi- 
tional difficulties arise studying biological fluids because 
of their complex rheological properties. 
In the cases of practical interest the flow is strongly 
viscous (creeping motions or Stokes flows) and a low 
Reynolds number regime can be assumed [3|. Since the 
creeping motion regime is a particular case of a lami- 
nar regime, it is possible to univocally define a time- 
independent pressure and velocity field. Therefore, for 
each position a well-defined drag force acts on a parti- 
cle immersed in the fluid flow. Ideally micro-flow sensors 
should be able to monitor the streamlines in real time 
and in the least invasive way. One common method to 
achieve this goal is to measure the drag force field act- 
ing on a probe particle, resorting to statistical criteria 
of analysis because of the intrinsic presence of Brownian 
motion. 

Recently optically trapped microscopic particles have 
been proposed as flow sensors 0, S S 0, S| ■ An opti- 
cal trap enables the confinement of micron sized objects: 
a high numerical aperture objective lens is used to tightly 
focus an optical beam and to produce a force sufficient 
to hold dielectric objects ranging in size from 100s of 
nanometers to 10s of microns [§]. In Q an oscillating op- 
tically trapped probe is used to map the two-dimensional 
flow past a microscopic wedge. In [B| and Q a stress mi- 
croviscometer is presented: it generates and measures 
microscopic fluid velocity fields, monitoring the probe 
particle displacement, which is directly converted into 
velocity field values, through digital video microscopy. A 
further improvement was achieved in [8] by using multi- 
ple holographic optical traps in order to parallelize the 
technique: an array of micro-probes can be simultane- 



ously trapped and used to map out the streamlines in a 
microfluidic device. All these techniques apply a Pho- 
tonic Force Microscope (PFM) approach (la Till. [liL Hi| 
to the flow measurement: the fluid velocity at the trap 
location is obtained by monitoring the probe displace- 
ments resulting from the balance between the trapping 
and drag forces. 

All the techniques described above present a major draw- 
back: they interpret the experimental results assuming 
that the flow can be described by a set of parallel stream- 
lines. When this is the case, the drag force acting on a 
probe around a specific point in the flow is well described 
by a constant value, and this is the case considered by 
the current optical trap methods. However, in the drag 
force-field there may be singular points as well. In these 
points the flow and, hence, the drag force are null, but 
not in its surrounding. The question, which arises natu- 
rally, is how to characterize this case. Under the assump- 
tion of a steady incompressible flow, a zero body force 
and low Reynolds numbers, the local fluid motion sat- 
isfies the quasi-static Stokes equations [3(: 7?V 2 v = Vp, 
V • v = 0, where v is the velocity, p is the pressure and 
rj the dynamic viscosity. Since the fluid is incompressible 
and there are no sources or sinks, there can be only two 
kinds of singular points [14j: (1) saddles (unstable) at 
the meeting point between two opposite flowing stream- 
lines or (2) centra (stable) in presence of a non-zero curl. 
As we show below, the characterization of the flow near 
these points can be achieved by studying the statistical 
properties of the Brownian motion of the probe. 
The knowledge of the flow-field near a singular point is 
of relevance for fundamental physical studies as well as 
for engineering applications. The mixing of fluids flow- 
ing through microchanncls is important for many chem- 
ical application; a reduction of the mixing length can be 
achieved by the generation of a transverse flow [l5| , case 
in which the formation of singular points is inevitable. In 
biological systems creeping motions take place in small 
blood or linfatic vessel or at the interface between tissues 
and prothesis or artificial organs [lfl]. In the presence 
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of slow flows, macrophage adhesion becomes more prob- 
able, increasing the possibilities of inflammation. Fur- 
thermore, the growth rate of thrombi, due to platelet 
aggregation, is also determined by the characteristics of 
the flow around it: for example, the presence of stable 
equilibrium poin ts, such as centra, helps their formation 
and growth 17] . Flow measurements at the microscale 
can help to diagnose pathologies and to guide the de- 
sign of bio-materials and nano-devices for diagnostic or 
therapeutic goals. 

In this letter we extend the PFM-approach to mi- 
crorheology in order to measure and characterize fluid 
fluxes in the proximity of singular points. The concept 
is to monitor the position of an optically trapped probe 
in order to locally characterize the drag force field as a 
generic function of the space coordinates up to the first 
order in its Taylor expansion around the probe position 
18]. This technique permits us to distinguish between a 
singular point in a complex flow and the absence of flow 
at all. Furthermore our approach allows one to deter- 
mine the stability of these singular points, which can be 
relevant for applications. 

In the following we will consider the drag force field 
produced by a bidimensional laminar regime. Since in 
microfluidic and lab-on-a-chip devices a planar geometry 
is generally assumed, this is the most useful case for ap- 
plications. Nevertheless, if needed, our approach can be 
generalized to the three-dimensional case. 

Experimentally, we generated these three kinds of fluid 
flow - namely a set of parallel streamlines, a saddle or 
a centrum - using solid spheres made of a birefringent 
material (Calcium Vaterite Crystals (CVC) spheres, ra- 
dius R — 1.5 ± 0.2/zm (H), which can be made spin due 
to the transfer of orbital angular momentum of light. 
They are all-optically controlled, i.e. their position can 
be controlled by an optical trap and their spinning state 
can be controlled through the polarization state of the 
light. The angular velocity of the CVC sphere is given 
by u> = AaAP/16ir 2 i]R 3 c, where Act is the change of the 
light spin momentum due to the scattering on the probe, 
P is the beam power, r\ is the medium viscosity, c is the 
speed of light, and A is the wavelength of the optical field. 
In the experimental realization up to four CVC spheres 
were optically trapped in water and put into rotation us- 
ing four steerable 1064 nm beams from a Nd:YAG laser 
with controllable polarization - to control the direction of 
the rotation - and power - to control the rotation rate. A 
probe polystyrene sphere (radius r = 500 nm) was held 
by an optical trap produced by a 632 nm linearly po- 
larized beam. The stiffness of this trap was adjustable 
through the power. The forward scattered light of this 
latter beam served to track the probe position using a po- 
sition sensor based on a quadrant-photodetector (QPD). 

Assuming no slip at the particle surface, the quasi- 
Stokes equation leads to the following solution for the 
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FIG. 1: (Color online) (a) Image of a trapped probe near a 
single spinning sphere. The bar is 1 \ivn. (b) 2D Probability 
Density Function (PDF) when the spinning sphere is at rest, 
(c-f) Total force fields - insets: drag force field contribution - 
and 2D PDF when the sphere is spinning (c-d) counterclock- 
wise and (e-f) clockwise. For each of the plots (b),(d), and (e) 
data acquired during 150 s with sampling rate 2000 Hz were 
analyzed. 



flow velocity near a single spinning sphere [H, Q 

R 3 



v = we, x x- 



(1) 



In each point near the rotating sphere the streamlines 
are perpendicular to the plane containing the z-axis (unit 
vector e z ) and the coordinate vector x. 

In Fig. 1 the results for probing such a drag force field 
near a single spinning sphere are presented. The CVC 
sphere is put into rotation at u> rs 38rad/s. The in- 
duced drag force can be measured through the shift of 
the equilibrium position of the probe (Fig. 1(b), 1(d), 
and 1(f)) and after calibrating the stiffness of the probe 
trap (k = 175 ± 15 fN//j,m). The corresponding force 
field can be reconstructed (Fig. 1(c) and 1(e)): the forces 
acting on the probe particle result 40 ± 5 fN (counter- 
clockwise rotation) and 37 ± 5 fN (clockwise rotation) . 
These results are in agreement with [j| @] . 

In a system with n spinning spheres, the linearity of 
the quasi-Stokes equations Q allows one to use the super- 
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FIG. 2: (Color online) (a) Force field and (b) auto- (red 
dashed line) and cross-correlation (blue solid line) functions 
when only the optical force is acting on a 0.5 \im radius probe 
in absence of drag force field. For plot (b) data acquired dur- 
ing 150 s with sampling rate 2000 Hz were analyzed. 



position of solutions (1) in order to obtain the total flow. 
Since the velocity in the equatorial plane is proportional 
to the reciprocal square of the distance from the sphere 
center, the flow velocity vanishes at a distance of a few 
sphere radii. Hence, in first approximation the hydrody- 
namic interaction between the spinning spheres located 
far enough from each other can be neglected. For n = 2 
or n — 4 we assume the spheres to be symmetrically dis- 
placed at the same distance p with respect to the origin 
x (Fig. 3(a) and 4(a)), so that the origin is a singular 
point, i.e. v(xo) = 0. The velocity field can be locally 
approximated as 



wi? 3 Jx = luR 3 



(2) 



where J is the Jacobian matrix of the total velocity field 
evaluated in the singular point and the Jacobian ma- 
trix of the velocity field generated by particle i evaluated 
in the same point. As expected, since the fluid is incom- 
pressible, Tr(3) is always null. This means that only 
(unstable) saddles - with n = 2, det(J) = — 4/p < - 
or (stable) centra - with n — 4, det(J) = 2/p > - can 
exist. 

fn Fig. 2 to 4 the experimental results for such cases 
are presented. We followed the data analysis proce- 
dure proposed in [13], to which we refer for the full 
details. The total force field acting on the probe is 
given by the sum of the drag force-field and the restoring 
force due to the harmonic trapping potential. The drag 
force field can be decomposed into a conservative and 
a rotational part, while the optical force field is purely 
conservative. The cross-correlation function (CCF) be- 
tween the x and y position of the probe, calculated as 
CCF{At) =< x{t)y{t + At) > - < y(t)x(t + At) >, docs 
not vanish only if the rotational part of the drag force 
field is not null. The rotational component of the drag 
force field can be obtained from the CCF 18, 191. If the 
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rotational component, and therefore the CCF, is null, the 



FIG. 3: (Color online) (a) Image of a trapped probe between 
a couple of spinning spheres. The bar is 1 \im. (b) 2D Prob- 
ability Density Function (PDF) when the spinning sphere is 
at rest, (c-f) Total force fields - insets: drag force field con- 
tribution - and 2D PDF when the sphere is spinning (c-d) 
counterclockwise and (e-f) clockwise. For each of the plots 
(b),(d), and (e) data acquired during 150 s with sampling rate 
2000 Hz were analyzed. 



total force field is purely conservative, and the conserva- 
tive component of the drag force field can be calculated 
by subtracting the optical potential from the total poten- 
tial. We underline that this can only be done if the total 
force field is conservative. Some results for an optically 
trapped probe in the absence of flow are presented in Fig. 
2: the optical force-field is harmonic (Fig. 2(a)) and it 
is purely conservative, as it is shown by the fact that its 
CCF function is null (Fig. 2(b)). 

In Fig. 3 the flow to be probed is generated putting 
into rotation (u> w 38 rad/s) two CVCs, symmetrically 
positioned with respect to the probe beads (Fig. 3(a)). 
The position of the probe is chosen so that there is no 
shift in the equilibrium position regardless of the rota- 
tion state of the spheres (Fig. 3(b), 3(d), and 3(f)). 
Since the experimental CCF is found to be null within 
the experimental error, we conclude that the rotational 
component of the drag force field must be negligible and 
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we can reconstruct its conservative part by subtract- 
ing the optical restoring force-field from the total force 
field (Fig. 3(c) and 3(e)). The optical trap produces a 
symmetric harmonic potential (k — 185 fN/ fim) (Fig. 
3(b)). With spinning CVCs the PDF becomes ellip- 
soidal (Fig. 3(d) and 3(f)). For a clockwise rotation 
(Fig. 3(c-d)), the main axes of the ellipses are oriented 
at 40° , the stiffness is kmin — 142 fN/ /im along the 
main axis and Umax = 261 fN/fxm along the secondary 
axis. For a counterclockwise rotation (Fig. 3(e-f)), the 
main axes of the ellipses are oriented at 42°, the stiff- 
ness is kmin = 165 ± 20 fN/fim along the main axis and 
kuAX = 211 ± 20 fN/fxm along the secondary axis. The 
difference are due to the fact that the rotation rate was 
observed to be different in the two directions. The cor- 
responding total force fields are represented in Fig. 3(c) 
and 3(c), while the drag force field contributions are de- 
picted in the insets: they constitute indeed saddle points. 
In Fig. 4 the characterization of a fluid flow near a singu- 
lar point of the second kind, a centrum, is presented. This 
stable singular point was generated putting into rotation 
(lo Ri 38 rad/s) four birefringent spheres symmetrically 
distributed with respect to the probe (Fig. 4(a)) in an 
optical trap with stiffness k = 78 ±5 fN/fxm. We can ob- 
serve that the PDF of the probe position do not depend 
on the rotation direction of the spinning spheres (Fig. 
4(b)); however, in this case the CCF is not null (Fig. 
4(d) and 4(f)), showing the presence of a rotational com- 
ponent of the drag force field. This component is such 
that it produces a torque on the probe bead, whose an- 
gular velocity is O, as jfl < x 2 + y 2 >= 5.7 fNfim for 
the beads rotating clockwise (Fig. 4(c-d)) and 6.4 fN/j,m 
for the case of the beads rotating counterclockwise (Fig. 
4(e-f)). The resulting total and drag force fields are rep- 
resented in Fig 4(c) and 4(e), while the drag force field 
contributions are depicted in the insets: they constitute, 
indeed, a centrum. 

In conclusion, we have demonstrated the use of an op- 
tically trapped probe for the characterization of singular 
points in microscopic flows. This technique delivers an 
important contribution towards a better understanding 
and optimization of microfluidic flows in the presence of 
singular points, providing both a deeper understanding 
of their physics and an experimental approach to their 
characterization. Furthermore our approach allows one 
to determine the stability of these singular points, which 
can be relevant for applications. 
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FIG. 4: (Color online) (a) Image of a trapped probe between 
four spinning spheres. The bar is 3 /im. (b) 2D probability 
density function when the spinning sphere is at rest, (c-f) 
Total force fields - insets: drag force field contribution - and 
auto- (red dashed line) and cross-correlation (blue solid line) 
functions when the sphere is spinning (c-d) counterclockwise 
and (e-f) clockwise. For each of the plots (b),(d), and (e) 
ten datasets acquired during 15 s with sampling rate 2000 Hz 
were analyzed and averaged. 
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